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Abstract. The Landau levels of cold atomic gases in non-Abelian gauge fields 
are analyzed. In particular we identify effects on the energy spectrum and density 
distribution which are purely due to the non-Abelian character of the fields. We 
investigate in detail non-Abelian generalizations of both the Landau and the symmetric 
gauge. Finally, we discuss how these non-Abelian Landau and symmetric gauges may 
be generated by means of realistically feasible lasers in a tripod scheme. 
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1. Introduction 

Gauge potentials are crucial for the understanding of fundamental forces between 
subatomic particles. A simple example of a gauge potential is provided by the vector 
potential in the theory of electromagnetism pQ. In this example the different vector 
components are scalars, and hence they commute with each other. If the vector 
components of the gauge field are not scalars, but instead N x N matrices, with N > 1, 
then it is in principle possible to have a situation where the different vector components 
do not commute. However, non-Abelian gauge fields are scarce in nature. Candidates so 
far have mainly been restricted to molecular systems [2] which are largely approachable 
only through spectroscopic means. Other systems are liquid crystals which show the 
required non-Abelian symmetries [3]. 

Experiments on cold quantum gases have reached an unprecedented degree of 
control, offering thus extraordinary possibilities for the analysis of the effects of gauge 
fields on atomic systems. A simple way of generating a gauge field in ultracold gases is 
obtained by rotating Bose-Einstein condensates with an angular frequency Qz (where 
we employ cylindrical coordinates {p, (p, z}). In the corresponding rotating frame 
the Hamiltonian describing the rotating system becomes one of a system subject to 
a symmetric Gauge field A = —mflpip, where m is the atomic mass [4J. Thus, a 
rotating condensate resembles a gas under the influence of a constant magnetic field 
B = mQ, and as a consequence many interesting phenomena, including e.g. Landau 
level physics and Quantum-Hall-like phenomena have been studied in rotating quantum 
gases [51 El Hi. 

Due to their internal structure, ultracold atoms offer as well the possibility of 
creating non-Abelian gauge fields. A surprising and astoundingly elegant derivation and 
description of the emergence of non-Abelian gauge potentials was presented by Wilczek 
and Zee [9]. It was shown by these authors that in the presence of a general adiabatic 
motion of a quantum system with degenerate states, gauge potentials will appear which 
are traditionally only encountered in high energy physics to describe the interactions 
between elementary particles. Ultracold atomic clouds are particularly promising 
candidates for realising such scenarios, since the access to physical parameters is, from 
an experimental point of view, unprecedented. Extending the ideas of Wilczek and 
Zee, it was recently proposed that properly tailored laser beams coupled to degenerate 
internal electronic states can be employed to induce Abelian as well as non-Abelian 
gauge fields in cold-atom experiments [TOj 121 [13] . Alternatively, such gauge potentials 
can be constructed in an optical lattice using laser assisted state sensitive tunneling 

!H nans El. 

With the implementation of these proposals, ultracold atoms would offer a unique 
testbed for the analysis of non-trivial effects on the properties of multicomponent cold 
atomic systems in the presence of non-Abelian gauge fields. To the best of our knowledge 
these effects have been scarcely studied in the literature [321 [13], [15l HHj. This paper 
is devoted to the analysis of non-Abelian effects on the spectral properties of ultracold 
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atomic systems. In particular, we show how purely non-Abelian effects lead to the 
eventual destruction of the Landau level structure, and may significantly modify the 
ground state density profile of ideal quantum gases. 

The structure of the paper is as follows. In Sec. [2]we study the generation of different 
forms of non-Abelian gauge fields, including non-Abelian constant fields, as well as the 
non-Abelian generalization of the Landau gauge. Sec. [3] is devoted to the analysis of 
constant non-Abelian gauge fields. Sec. 0] discusses the non-Abelian Landau gauge, 
and in particular the destruction of the Landau level structure and the corresponding 
modified de Haas- van Alphen effect. In Sec. [5] we discuss the non-Abelian symmetric 
gauge. Finally in Sec. [6] we conclude and discuss some future promising directions. 



2. Laser-induced non-Abelian Gauge fields 

In this section we discuss the generation of (possibly non-Abelian) gauge fields as those 
discussed in the following sections. There are at least two alternatives for the creation 
of non-Abelian gauge fields. One consists in employing two-component atoms in state- 
dependent optical lattices in the presence of appropriate laser arrangements [HI [15]. 
A second possibility, which we shall explore in this paper, was recently proposed in 
Ref. [12]. 

In this second alternative, a non-Abelian gauge potential is constructed for 
atoms with a tripod electronic structure [TjJJ [20]. Three lasers with appropriate 
polarizations couple the excited electronic state |0) and the groundstates \j = 1,2,3), 
with corresponding Rabi frequencies f^-(r) which parametrize in the form: Qi = 
Q sin 6 cos</>e , f2 2 = sin# sin <p e tS2 , Q% = Q cos#e iS3 . For a fixed position r the 
Hamiltonian describing the laser-atom interaction may be diagonalized to give a set of 
dressed states. Under appropriate conditions two dressed states, so called dark states 
become decoupled from the other states: 

|.Di) = sin0e i531 |l) - cos0e i532 |2), (1) 
\D 2 ) = cosflcos0e i531 |l) + cos#sin0e i532 |2) -sin0|3), (2) 

with Sij = Si — Sj. The dark states have zero eigenvalues and are separated by the energy 
HQ from the remaining eigenstates. The state of the atom can therefore be expanded in 
the dark state basis as |$) = \E f 1 (r) \Di(r)) + ^(j) l-^tj))- The two-component spinor 
^ = ^2} T obeys a spinor Schrodinger equation of the form 



k 9 T 



—{-ihv - Ay +v + $ 

2m 



*• (3) 
In this equation we observe the appearance of a 2 x 2 vector potential of the form 
An = ft (cos 2 (f)VS 2 3 + sin 2 0VSi 3 ) , 

A12 = ft COS 9 (isin(20)VS 12 -iV0) , (4) 



2 v r, ~- -rj 

A22 = h cos 6 ( cos </>VSi 3 + sin 2 4>VS- 
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The systems also presents a scalar potential of the form $jj = • kj, where 

Ki = sin 9 Q sin(20) V5i2 + iV^j , (5) 

k 2 = - sin(2fl) (cos 2 4>VS 13 + sin 2 0VS 23 ) - iV9. (6) 

Finally, if the original states \ j = 1,2,3) experience an external potential Uj(r), then 

V n = U 2 cos 2 <p + C/i sin 2 0, (7) 

Vi 2 = ^^ cosflsin(20), (8) 

V 22 = (Ux cos 2 + U 2 sin 2 0) cos 2 9 + U 3 sin 2 0. (9) 

Recent advances in shaping both the phase and the intensity of light beams make 
it possible to achieve a remarkable versatility in controlling the gauge fields [21, 22J, 
provided the corresponding light fields obey Maxwell's equations. 

In the following we shall assume that the atoms are strongly trapped in the z- 
direction, hence they are confined to the xy plane. Given two orthogonal vectors £ and 
ff on the xy plane, we shall be interested in non-Abelian situations, in which A^ = A ■ £ 
and A v = A ■ ff, fulfill [A^, A v ] 7^ 0. This condition demands (u x V<f)) z ^ and/or 
(u x VS 12 ) Z ^ 0, and/or (VS 12 x V0) z ^ 0, with u = (cos 2 - cos 2 9 sin 2 9) VS 23 + 
(sin 2 - cos 2 9 cos 2 (j))VS 13 . 

2.1. Constant intensities 

We will consider first homogeneous intensity profiles, i.e. both and 9 are now space 
independent. We choose the particular case with = 9 = tt/4. For constant the 
non-Abelian character demands VS23 x VS\ 3 7^ 0. A simple laser arrangement fulfilling 
this condition is Sj 3 = ajx + j3jy, where Oj,ft are constants such that a 2 ft 7^ aift. 
The corresponding x and y components of the vector potential are of the form 

A x = ~(ai + a 2 )(3i + a z ) + _ "2) £3, (10) 

K = + ft)(31 + a z ) + ^=(ft - (3 2 )a x . (11) 
On the other hand, by choosing Vj(r) = AEj + U(f), with AEi = — (h 2 /16m)[(a 2 — 



a 



I) + (/3f - ft 2 )] = -AE 2 , and AE 3 = -(ft 2 /16m)[(a 2 + a\) + + ft 2 )], one can prove 
that (up to an irrelevant constant) V+(f) = U(f) with U(r) a common trapping potential 
for all components. 

A gauge transformation eliminates the terms proportional to the identity matrix in 
A x and A y . Let hn y = (ft - /3 2 )/2y/2, hq y = (ft + ft)/8, Hk x = {a x - a 2 )/2\/2, and 
hq y = («i + a 2 )/8. A rotation a x — > cosr]a x + smr]d- z , a x — > — smna x + cosna z , 
with tan2r/ = K y /q y , provides A y = hq y d z , with q y = cos20g y + sin20K y , and 
A x = hk x a x + hq x a z , with k x = (cos 20/^ — sin 20^), and q x = (cos20g x + sin20fi; z ). 
Hence, we recover exactly the same form which was discussed in Sec. |3j 
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2.2. Landau-like gauge 

In this subsection we shall consider the case S13 = S23 = S. In that case the non-Abelian 
character demands (VS x V0) 2 7^ 0. We will choose the phase S = nx, and <fi = qy, 
which gives a non-Abelian gauge potential unless k = or q = 0. In addition we take 
cos 6* = x/R c , where R 2 , = x 2 + (z — z c ) 2 , such that for the relevant x-range, \x\ « z c 
is fulfilled. As a consequence, and up to first order in (x/z c ) we obtain: 

A ~ hn(l + d z )x + B xa y y. (12) 

where Bq = — . Note that although x <C z c , B can actually have large values. In 
addition, and again up to first order in (x/z c ), we obtain V + (ft = U(r), if V\{r) = 
= h 2 q 2 /2m + U(r), and V${r) = h 2 /2mz 2 . Using a simple gauge transformation 
\& — > expinxty to eliminate the identity matrix term in A x , and applying a unitary spin 
transformation WAlI, with U = (a z + a y ) / y/2, we obtain A ~ fiKa y x + B xa z y, which 
is indeed exactly the same Landau-like gauge that we employ in Sec. IU A simple laser 
arrangement which would lead to this particular gauge is provided by 

fli = tt cos qye iK{x+y+z)/2 , (13) 
Q 2 = Qsinqye l ^ x+y+z)/2 , (14) 
Q 3 = n^e^-^/ 2 , (15) 

Z c 

where we assume the illuminated atoms are confined to a region for which \x\ <^ z c 
holds. 



3. Constant non-Abelian gauge 

Let us consider a constant matrix gauge of the form A = (A x , A y , 0). We have already 
shown that these fields can be generated in a tripod scheme using a simple laser 
arrangement. Then, the Hamiltonian of the 2D system becomes: 



1 

H = — 

2m 



(P x + A X f + fa + A y f . (16) 



In the Abelian case [A^Ay] = 0. We can therefore choose a common eigenbasis for 
both matrices: A x /h = diag{gi x ., g 2x }, A y /h = diagjg^, q 2y }, As a consequence, we 

recover two independently displaced quadratic spectra Ej(k) = ^ (k + qj , where 

Qj = [Qjx, Qjy)- 

In the non-Abelian case, on the other hand, we cannot simultaneously diagonalize 
both matrices, and as a consequence the spectrum becomes distorted. Let us consider 
a simple, but representative, case, namely A x = q x a x , A y = q y a z . Employing the 
Fourier-like transformation 

Vfc, y) = J2 e ik ^ (^p) e ik ^$(k x , k y ) (17) 

kx,ky v V2 y 
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with k XjV = 27rn X) y/L, we may transform the time-independent Schrodinger equation 

Eijj(x,y) = Hip(x,y) into 

2mE -* -* -* 

—^-4>{k x , k y ) = [k 2 x + q 2 x + {k y + q y ) 2 ] 4>(k x , k y ) + 2q x k x 4>(k x , -k y ). (18) 

Diagonalizing the system of equations for 4>(k x , ±k y ), we obtain two eigenenergies 



2mE. 



± _ 7.2 , „2 , j 2 , 2 



K + Qi + K + ly^^^ + Vv ( 19 ) 

Note that in the Abelian situation q x = (or q y = 0), and, as expected there is no 
coupling between momenta in different directions. However, due to the non-Abelian 
character, even for a constant gauge there is a non trivial coupling between the different 
directions. 



4. Landau-like Non-Abelian gauge 

4-1. Periodic boundary conditions 

We consider in the following a matrix generalization of the Landau gauge, namely 
A = (hnM x , B M y x, 0) (the usual Landau gauge is of the form (0,B x, 0)). We will 
assume that the matrices M x and M y are constant. Then the Hamiltonian of the 2D 
system becomes: 

2m 



(p x + HkM x ) 2 + (p y + B Myx) 2 . (20) 



We first discuss the typical text book situation, in which the particles (which are 
assumed to be confined on the xy-plane) are considered as confined in a 2D box of side 
L with periodic boundary conditions (i.e. a toroidal configuration). We are particularly 
interested in how the non-Abelian character of the fields destroys the usual Landau- 
level structure of the energy eigenstates. In the following subsection we shall discuss a 
slightly different scenario closer to actual experimental conditions. 

As in Sec. |3l if [M x , M y \ = 0, one can find a common eigenbasis {e*i, e*2}, such that 
in this basis M x = diag{7i,72}, and M y = diag{Ai, A2}, and hence the Hamiltonian is 
also diagonal in this basis. Since we assume periodic boundary conditions we can thus 
consider wavefunctions of the form 



^•(f) = e'^SK (21) 



such that 



E Vj (q) 



P 2 . 1 22 
— muj.q 

2m 2 jH 



vAq). (22) 



where q = x + ^ , V = —ihd/dq, and tUj = B \\j\/m is the cyclotron frequency for 
the state j. Hence, for the Abelian case we obtain two different sets of Landau levels 
with energies Ej(n) = t\Wj{n + 1/2), and degeneracies gj = B XjL 2 /2iih. Note that 
if I Ai| = I A2 1 5 as it is the case for M y = a z , then the two sets of Landau levels are 
degenerate. 
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Let us now discuss what happens if on the contrary [M x , M y \ 7^ 0. We work (without 



lack of generality) in the basis in which M y = a z . Note that the Ansatz 



(23) 



also fulfills periodic boundary conditions, 
equation to obtain 



We insert this Ansatz in the eigenvalue 



Eu(n y , x) 



+ 



2m 2m 



d- z [fl 2 ,a z 



2im. 



- + T X 



u(n y , x) 



u(n y , x) — u(—n y , x) 
4m 



(24) 



where IL = p x + hnM x . For the Abelian case, [M X) a z \ = 0, the last term vanishes, and 
we get the same equation as previously. However, for the non-Abelian case, the last 



term introduces a coupling between the modes with n y and 



-n 



y 



and hence there is 



an explicit dependence on n y . As a consequence of that, the degeneracy of the Landau 
levels is lifted. 

For the particular case of M x = a y , we get the following set of coupled equations 
(e = E -h 2 K 2 /2m): 



eu(n y , x) 



Pi 



Bl 



2m 2m 



x + 



2ixtm % 



u(n y , x) 



+ —p x d- y u(-n y ,x) 
m 



(25) 



eu(—n y , x) = 



f) 2 B 2 
Px _j_ £^o_ 

2m 2m 



x 



27ihn l 
B n L 



u(—n y , x) 



H PxO-yU{n y ,x). 

m 



(26) 



The coupling prevents the re-absorption of n y in the definition of a new q variable, as it 
was done in the Abelian case, and hence the spectrum explicitely depends on n y . Note 
that we are imposing periodic boundary conditions, and hence a; is in a ring of perimeter 
L. In this sense, ±L/2 are the same point, and this must be taken into account when 
considering the harmonic oscillator potential in each equation. 

Note that the previous equations involve the coupling of harmonic oscillator 
wavefunctions centered in ±x c (\n y \), with x c (\n y \) = 27ih\n y \/ B L. Hence, the smaller 
the overlapping between coupled wavefunctions, i.e. the large x c , the smaller the 
coupling, and as a consequence only sufficiently small values of n y will be affected by the 
non-Abelian coupling. This point becomes clear after performing first order perturbation 
theory assuming a small coupling k. A straightforward calculation shows that the lowest 
Landau levels, which correspond to the lowest eigenvalues of each harmonic oscillator, 
experience a maximal energy shift 

- tl ( /c g^_ e -nS/Ar^ > (27) 



An, 
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where l\ = h/muo c is the magnetic length, and An y = y/g/2ir, with g the degeneracy 
of the unperturbed Landau levels. Note that for n y = the first correction should 
be quadratic in k, whereas for n y ^ it should be linear. Clearly, the relative 
importance of the non-Abelian corrections should decrease as \j ^fg. In particular, 
the maximal energy shift (AE) averaged over the different n y can be approximated as 
{AE)/hw c ~ ( K l c )/^g~. 

We have solved numerically for the eigenvalues of Eqs. (1251) and fl26|) imposing 
periodic boundary conditions, for different values of g which controls the strength of 
the magnetic field applied, and kI c which provides the strength of the non-Abelian 
corrections. The value of L/l c = \J1-ng is chosen in all simulations. Fig. [1] shows 
the behavior of the lowest eigenvalue as a function of n y for g = 128 and kI c = 
0,0.2,0.4,0.6,0.8, 1.0 (from the uppermost to the lowermost curve). The figure follows 
approximately the perturbative result. For n y = & higher order contribution appears, 
but note that a quadratic law follows for small k follows, and not a linear one, as in the 
case for n y ^ 0. As expected from the previous calculations only values of n y up to the 
order of ^fg contribute significantly to the shift of the lowest Landau level. 




fly 



Figure 1. Lowest eigenvalue e/Hui c as a function of n y , for g — 128 and kI c = 
0, 0.2, 0.4, 0.6, 0.8, 1.0 (from the uppermost to the lowermost curve). 

Figs. 2 show the behavior of the Landau levels for g = 128, and kI c = (a), 0.6 (b). 
The figures are presented as histograms in intervals of 0.05hu c , in order to reveal more 
clearly the destruction of the Landau levels. Note that the gaps (of energy fku c ) between 
the Landau levels are filled, and the peaks in the density of states are progressively 
reduced. For sufficiently large k the Landau level structure therefore disappears. 

4-2. Absorbing boundary conditions 

In the previous section we discussed how the non-Abelian character of the gauge field 
significantly modifies the text-book Landau level structure. In the following, we consider 
a slightly different physical scenario which is closer to the actual experimental conditions 
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Figure 2. Landau level structure for periodic boundary conditions, g = 128, and 
kI c = (top) and nl c = 0.6 (bottom). We employ (see text) M x — a y and M y = a z . 



discussed in Sec. |2j The particular procedure devised for the generation of the non- 
Abelian Landau gauge demands that the x coordinate cannot be considered as periodic. 
We take the same box configuration as for the previous subsection, but assume absorbing 
boundary conditions in the x direction, while keeping for simplicity periodic boundary 
conditions in the y-direction. We consider exactly the same gauge discussed in the 
previous subsection. The spectrum is provided by Eqs. (125]) and (1261) but imposing 
absorbing boundary conditions. Figs. 3 show the lowest Landau levels for the same 
cases discussed in Figs. 2. 

Even for the Abelian case the Landau level structure is of course affected by the 
absorbing boundary conditions. In the Abelian discussed in the previous section, 

the problem reduces to two decoupled equations for harmonic oscillators centered at 
±x c (\n y \). Clearly when x c approaches L the levels of the resulting potential become 
largely distorted, leading to a significant modification of the Landau level structure when 
n y approaches g. This reduces the effective degeneracy of the lowest Landau levels to 
values smaller than g. The effective degeneracy, as shown in the figures, becomes smaller 
for higher Landau levels. The non-Abelian effect leads, as in the previous subsection, 
to the eventual destruction of the Landau level structure. 




Figure 3. Same cases considered in Fig. [5] but with absorbing boundary conditions 
(see text). 
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4-3. Modified de Haas-van Alphen effect 

The destruction of the Landau level structure has experimentally relevant consequences 
for the behavior of cold atomic gases. As an example we can consider the case of 
an ideal two-component Fermi gas under the previously mentioned non-Abelian gauge 
potential (we consider a temperature T <C Tp, where Tp is the Fermi temperature). 
Equivalently to the well-known de Haas- van Alphen effect [23] , we may study the energy 
per particle, E = E/N, of the Fermi g function of the applied magnetic field 

B , or equivalently of g. This energy may be monitored by measuring the released 
energy in time-of-flight experiments. For k = (Abelian case) d 2 E(B )/dBQ presents 
a typical configuration of plateaux, due to the degeneracy of the Landau levels. The 
destruction of the Landau level structure significantly distorts this picture, rounding-off 
or eventually destroying this plateaux configuration (see Figs. Hj). 
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Figure 4. Value of —d 2 E/dB^ as a function of the applied magnetic field Bq, for the 
same case discuss in Fig. [U and k — (filled circles) and k = 0.6 (hollow circles). 



5. Symmetric gauge 



In this section we consider an ideal cold atomic sample in an isotropic harmonic trap 
of frequency u, in the presence of a non-Abelian generalization of the symmetric gauge 
of the form A = A p p + pA v 0. Although the tripod scheme is not suitable for the 
experimental realization of this gauge, we include the analysis of this gauge field for 
completeness of our discussion. Other ways of generating non-Abelian gauge fields, as 
lattice techniques [T5] should be employed in this case. In the following we consider 
A p = HnUp, A v = BqU^, where U p>tp are linear combinations of {1, a x , a y , a z }. 
The corresponding time-independent Schrodinger equation is of the form 



Eip 



1 

2m 



-ihV + A 



mu 



ip H — p tp. 



(28) 



Performing the gauge transformation ip 
transforms into 

1 

2m 



exp[— iA p p/h](j), the Schrodinger equation 



Eip 



-ihV + <pC v (p)p 



Tp + 



mu 



-P 2 *P, 



(29) 
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where 







where we reduce the equations to a dimensionless form by employing oscillator units 



cyclotron frequency. 

As mentioned above the non-Abelian character of the gauge field induces an 
additional p-dependent potential. It severely distorts the standard Fock-Darwin 
spectrum which is expected for the Landau-level structure in the presence of a symmetric 
gauge and a harmonic potential, as it is shown in Fig. [5j An inspection of the level 
structure shows that not only the eigenenergies are modified, but also the ordering 
of the different eigenstates becomes distorted as a consequence of the non-Abelian 
potential. As a consequence of this extra p-dependent potential, an ideal Fermi gas 
at zero temperature shows a significantly distorted density profile in the presence of the 
non-Abelian gauge field, as shown in Fig. [6j 

6. Conclusions 

In this paper we have analyzed the physics of ultracold gases in the presence of a non- 
Abelian gauge field. We have first studied how different types of non-Abelian fields may 
be created by means of relatively simple laser arrangements with atoms described by 
an electronic tripod level scheme, including a non-Abelian generalization of the Landau 
gauge. In a second part we have considered the non-trivial effects that the non-Abelian 
character has on the eigenlevel structure of the cold atomic system. In particular we 
have shown that exclusively due to the non-Abelian character of the field, the usual 
Landau level structure is severely distorted, and even eventually destroyed. We have 
shown that this effect may be observable in an equivalent experiment to the well-known 
de Haas- van Alphen effect. The distortion of the Landau levels leads to a significant 
modification of the usual plateaux-like signal characteristic for the de Haas- van Alphen 
effect. Finally, we have completed our analysis of a non-Abelian version of the symmetric 
gauge. We have shown that the Fock-Darwin spectrum is significantly distorted in the 
presence of non-Abelian fields, due to the presence of an extra potential, which is a 
purely non-Abelian effect. 




. In the previous equation 
Bq/ttl is the corresponding 
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Figure 5. Fock-Darwin spectrum E/Hlu for the Abelian (top) and non-Abelian 
(bottom) cases discussed in the text as a function of bo = uj c /llJ. 



30- 




Figure 6. Comparison between the density profile for an ideal Fermi gas occupying 
up to 56 eigenlevels at zero temperature for the Abelian (solid) and non-Abelian cases 
discussed in the text with different values of k = 1 (dashed), and n = 5 (dotted). 
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